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Abstract. In this work, we present a stability criteria for the solitary wave 
solutions to a BBM system that contains coupled nonlinear terms. Using the 
idea by Bona, Chen and Karakashian [5] and exploiting the accurate point 
spectrum information of the associated Schrodinger operator, we improve the 
stability results previously got by Pereira [15]. 


1. Introduction 

This work is essentially motivated by two papers, [5] and [15]. In [5], Bona, Chen 
and Karakashian consider a class of BBM-type systems in form of 


( 1 ) 


u t +u x - u xxt + P(u, v) x = 0, 
V t + V x - v xxt + Q(u, v) x = 0, 


where 

P{u, v ) = Au 2 + Buv + Cv 2 , Q(u, v) = Du 2 + Euv + Fv 2 

with constants A,B,--- ,F £ R. Namely, P and Q are homogeneous, quadratic 
polynomials. The other paper [15] by Pereira considers BBM system 

(2) Ut + co U x — AU xxt + (V H(u)) x = 0, 

where U = U(x, t) is an 7?. 2 -valued function, Cq is a non-negative parameter, A is a 
2x2 real positive definite matrix, and VH is the gradient of a homogeneous function 
H : 1Z 2 — > 1Z with proper regularity. The nonlinear terms in (2) can be seen as 
generalizations of that in (1), since the latter were essentially forced into a structure 
of the same form as the former, due to certain positiveness requirement. However, 
the result in [5] turns out to be stronger than that in [15]. The reason is that, in 
a proper space, [5] takes into account the accurate point spectrum information of 
the Schrodinger operator £ 0 = — -js — where </> is a typical potential function. 
Here we present a stability result which generalizes [5] and improves [15]. 
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1.1. Main result. We consider a BBM system in form of 

(3) U t + U x - U xxt + ( VH(U)) X = 0, 

where U = U(x,t) is an TZ 2 valued function and WH = ( H Ul H v Y is the gradient 
of a C 3 homogeneous function H : TZ 2 —> 1Z , the superscript * represents transpose. 
The equations are the same as in (2), except that we have simplified the system by 
letting co = 1 and A = 12x2 for sake of clarity. The techniques used in this work 
apply to cases with more general co or A. In this note we focus on the stability 
results. We refer to [9, 14] for global wcll-posedness of (3), or [15, 9] for instability 
results. 

The functional space is L 2 (JZ ), in which the scalar operator C defined below has 
domain H 2 (1Z). While in the coupled system, the natural product spaces L 2 (JZ) x 
L 2 (1Z) and H 2 (1Z) x H 2 (1Z) are used. The operator C sends H 1 ( TZ) to H~ 1 (TZ); 
< Cu,u > denotes the the pairing of Cu and u on space H~ 1 (TZ) and H 1 { TZ). We 
shall not further dig into the issues on space setting, instead we simply assume all 
the functions here have enough regularity for all the operations under consideration. 
However, precise functional settings are crucial in both the wcll-posedness and the 
stability proofs, for instance, in case of the space decomposition [6]. 

We assume here the homogeneous function H(-,-) is of order p + 2 with integer 
p > 1. Suppose <f> = (</>, fMf>y is a proportional solitary wave solution of (3), where 
p G TZ is the proportional coefficient. Assume H u (l,p) > 0 and let 

a a _ _1_ ( H uu (l,p) H uv { 1, p) \ 

{p T l)Hu(l, p) \ H uv (l , p) H vv (l, p) J 

Our main result is 

Theorem 1.1. Assume det(Al) < ^j-j-. We have 

(1) for p < 4, $ is stable; 

(2) if p > 4, then there exists a u> p > 1 such that is stable for uj > uj p and 
unstable for 1 < ui < u> p . 

Remark 1.2. Secret-telling moments. 

(1) In this note, we only consider the proportional solitary wave solution. Nev¬ 
ertheless, as the study of nonlinear Schrodinger system suggests [19], non¬ 
proportional solitary wave solutions do exist and are of great importance in 
real applications. 

(2) Given H, one can easily compute det(A4) and determine the stability of a 
solitary wave. This idea was first brought up in [5]. The accurate infor¬ 
mation on eigenvalues is the key to discover the condition under which the 
stability is achieved. 

1.2. Strategy. Since Benjamin [3] set up the rigid mathematical framework, with 
refinements, generalization and elaboration along the way [4, 16, 17, 18, 8, 2, 6, 13] 
etc., stability study in this direction has developed into a sophisticate but “routine” 
process. On the one hand, it has close relations with areas such as scattering/inverse 
scattering theory, spectrum theory, hamiltonian systems, etc., see [12]; on the other 
hand, the classical method more or less follows a routine process in proving the 
stability of a solitary wave solution, <f> for instance: 

• formulate the proper metric 11 • 11 to measure the shape deviation; 

• capture a pair of conserved quantities to form a Lyapunov functional L , 
which is invariant too; 
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• for any solution u with perturbation v = u — </>, investigate L(u ) — £(</>) and 
show there exist positive constants C \, C 2 such that for all time, 

(4) CilMI 2 + o(|M| 2 ) < L(u) m < C 2 \\v\\ 2 + o(||u|| 2 ). 

Then the second inequality at time t = 0 implies that L(u) — L(<f> ) is bounded 
by initial deviation ||u( 0)|| 2 = ||ix( 0 ) — (f || 2 ), at least when v is small; the first 
inequality at all later time implies all deviations ||v || 2 are bounded by L(u) — £(</>). 
The fact that L(u) — L(cf>) is invariant gives the stability: the initial deviation (shape 
difference) bounds the future deviation, at least when the initial deviation is small. 

This clear picture does not show up without subtle issues. The global well- 
posedness needs to be ready; the key player-a solitary wave solution has to exist; 
to maintain an effective norm, the infimum of the deviation has to happen within a 
finite translation [4]; how to transfer a constrained result to a general one, etc. But 
the most delicate part is the establishment of inequality (4). 

In estimating L(u) — a self-adjoint operator £ shows up. Roughly, we 

have L(u ) — L(<f>) =< Cv,v > +o(||u|| 2 ), where < •,• > denotes the adjoint pair 
on H~ 1 (TV) x H 1 (TV). The upper bound, L{u) — L(<j)) < C 2 IHI 2 + o(||u|| 2 ), is 
easy to show, since £ is bounded operator from H 1 ^) to H~ 1 ( 1Z). The lower 
bound requires certain amount of “positiveness” of £. The operator is generally 
not positive definite, but almost so: the residual spectrum is missing since the 
operator is self-adjoint; it has strictly positive continuous spectrum; it has a simple 
zero eigenvalue and a unique, simple, negative eigenvalue. Thus, the space under 
investigation is able to be decomposed into a direct product of certain subspaces 
and the major part of the positiveness is kept, so is the positiveness of the lower 
order term ||u|| 2 . 

This classical method calls for a pack of techniques that attach to the analysis of 
such an operator. To avoid these difficulties, the concentration-compactness method 
has been developed by Cazenave and Lions from the calculus of variation viewpoint, 
see [7, 11]. 

We shall embrace the tradition and follow the process mentioned above. But 
as has been remarked before, we shall omit some “critical but routine” steps and 
focus on the those parts which reflect the improvements we provide: the accurate 
spectrum analysis on the Schrodinger operator; the calculation of a critical quantity 
Some examples are given in the end. 

1.3. Preliminaries. 

1.3.1. Point spectrum of Schrodinger operator. The result in this work relies on 

the precise information on the point spectrum of a specific Schrodinger operator, 

,2 

£0 = —A + V(x) in one dimensional space, where the Laplacian A = ^3 and V 
is a potential function. In the special case where the potential function Vix) = 
—asecli ( 2 ;), a > 0, the eigenvalue problem can be solved explicitly; see Landau 
etal. [10], page 73-74. Here for sake of completeness, we reproduce some details in 
this aspect. 

On space L 2 (1Z ), this <j>{x) make £0 a self-adjoint operator. Let a denote the 
full spectrum and let o p , a c , oy be the point spectrum, the continuous spectrum and 
residual spectrum, respectively. For self-adjoint operator the residual spectrum is 
missing, we have cr(£ 0 ) = cr p (£o) U cr c (£o). Moreover, 

Lemma 1.3. [Pava [12]] we have 

(1) has spectrum a = a c = [ 0 ,oo). 
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(2) — fpps + u) has spectrum a = a c = [uj, oo). 

(3) Cq = —^2 + u) — asech 2 (x), being self-adjoint on L 2 (JV), has only real 
spectrum, which consists of the continuous part [u;, oo), together with a finite 
number of discrete eigenvalues in (—oo,w). 

To find the eigenvalues of C 0 , we shall solve 


(5) 


d 2 2 

—~ a sech ( x)ip(x ) 




Make variable transformation f = tanh(x), hence £ G (—1,1) for x € (— 00 , 00 ), 
we have the the associated Legendre equation for transformed function fj(f) = 

V>(tanh _1 (£)) : 


( 6 ) 


d 

df 




A 

1 — £ 2 


M) = 0. 


Note that we work in space L 2 (1Z), which require ip{x) —> 0, as x —> ± 00 , or 


■ 0 (£) “^ 0 , as £ —> ± 1 . 


Letting 


(7) e = Vx, s(s + 1) = a, s > 0, (hence s = ^ ^ — 1 + \/l + 4aj) 

and putting the equation in form of 


( 8 ) 


d 

df 


U-C 2 )-^ i>(0+ s(s + i) 

df J L 


1-C 2 


^(0 = o , 


we can solve the eigenvalue problem ( 8 ) completely: 

X n e n , e n s |Y| ~b ^ 1,2, • , [”s~|, 

(9) MO = p/»(0, 


where |~s"| is the smallest integer greater than or equal to s, and is are the 

associated Legendre functions of the first kind , meaning the functions are defined 
within the unit circle in the complex plane. By analytical continuation, Legendre 
functions of second and third kind can be defined outside of the unit circle. But 
they do not interest us here. 

When s is a positive integer, e = 1,2, ■ • ■ , s and Pf becomes the Legendre poly¬ 
nomial which has an explicit form, for example when s = 3 


Pi(0 = 3(1 — £ 2 ) 1/2 (15£ 2 — 3), 

Ps(0 =15(1-a?, 

P|(0 = 15(1 — £ 2 ) 3/2 . 
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The solutions for the corresponding eigenvalue problem (5) are 


Ai 

= - 1 , 



= P 3 1 (tanh(x)) 

= 3/2(1 — 5tanh 2 (:r))sech(a;); 

A 2 

= — 4, 


t/) 2 (x) 

= P 3 (tanh(a;)) 

= 15tanh(x)sech 2 (;r); 

A 3 

= -9, 


i>3{x) 

= P 3 (tanh(x)) 

= — 15sech 3 (x). 

not an 

integer, the expressions of P/(£) can be put in form of hyper- 


geometric function (page 74, [10]) 


( 10 ) 

whe 

( 11 ) F[-m,P,'y,z] = 


P e s n (?) — (1 — ? ) 2 F[e n — s, e n + s + 1 , e n + 1 , — (1 — ?)], 


where the hyper-geometric function F is defined as 

Z 1 -T(l - z)T+™-/3 d” 


7+m —1 


(1 -zf-' 1 ]. 


7(7 + 1) • • • (7 + m — 1) dz r ‘ 

Again here what concerns us is the information of the concrete eigenvalues, rather 
than the eigenfunctions. 

Another result will be used later. 

Lemma 1.4 (Theorem B.61,[12]). Suppose that <f> 6 L 2 (TV) satisfies the differential 
equation 

-0" + ccj) - f{ 4 >) = 0, 

with c > 0 and (/>' having exactly one unique zero. Then the differential operator 

Lo = -^ + [c ~ f ' m 

defined in L 2 (JZ) has exactly one simple negative eigenvalue Xq; the eigenvalue 0 
is simple with associated eigenfunction <f>; and there exists S > 0 such that every 
X £ ct(Lq) — {Ao, 0} satisfies A > S. 

1.3.2. Scalar equation. In this part, by comparing with the scalar equation, we 
single out the difficulties brought up by the coupling mechanism and focus on these 
new issues. 

Consider a scalar BBM equation 
(12) u t + u x — u xxt + Hu p u x = 0, 

where h > 0 is a constant. There exist infinitely many conserved quantities, two of 
which play an important role in our analysis: 

1 /' ■ 


^ o ( m ) = 7 : I u~ 

z Jn 

©o(w) = -~2 J n u2 + 


iidx. 


2 H 


-u p+2 dx. 


2 Jn (p + 2 )(p+l) 

We assume the solitary wave solution exists in form of <j) u (x — cot). Then let 
z = x — tot, (f>{z) = 4> u {x — cot), substitute into ( 12 ), we have 

-cocfi + (f' + ucfi" + {-^-cf p+1 y = 0 . 

P + 1 
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Requiring (j>{z) decreases to zero at infinity, we have equation associated with 
( 12 ) 

(14) — (w — 1 )(f> + uj(f>"-\ - -(f) p+1 = 0. 

p + 1 

When u> > 1, the solution to (14), which is 


</> 0) = 


(p +2)^-1) 2 p 

v 2 


2 h 


u> — 1 


w 


*) 


generates the solitary wave solution 4> u (x — cot) to (12). 
Taking derivative of (14), we have 


1 

V 


— (w — 1 )ft + toft" + heft ft = 0 , 


from which we conclude the operator 


C 


sclr — 



(w — l) — hfP 


has eigenfunction ft corresponding to eigenvalue 0. 

In view of w > 1 and (f> being in shape of [sech 2 (-)]p, the above is in fact the 
consequence of Lemma 1.4. Furthermore, we conclude 

1. Csdr has exactly one simple negative eigenvalue Ao; 

2. the eigenvalue 0 of C sc i r is simple; 

3. and there exists <5 > 0 such that every A £ a(C sc i r ) — {Ao, 0} satisfies A > S. 
Let 

d(w) = 0 O 


The classical result is, the main inequality (4) will be established if 


d"(w) > 0 


holds, however, under the constraint 

= f2o(^w) = Const. 

Then an extra argument (omitted in this note) shall bridge this constrained case to 
the general case and Theorem 1.1 would then be completed. 


Remark 1.5. For a general solitary wave solution, there is a requirement on the 
existence of an interval of uj, on which d(ui) has certain smoothness. It is not a 
problem for the special solitary solution here. 

Four the coupled system (3), we shall walk through the same procedure, however, 
with conserved quantities and operator in different forms. And the main difficulty 
comes from the spectrum analysis of the new operator. 


1.3.3. Coupled system. Consider (3), 

Ut + U x - U xxt + C VH(U)) X = 0, 
We have conserved quantities 

V(U) = \ j 11*11 + UlU x dx, 

z Jn 

0(17) = —i [ U t U + 2H(U)dx. 
2 Jn 
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We consider the proportional solitary wave solution U(x,t) = — wt) = 

(4>ui{x — u>t),if) u (x — uit)Y of system (3), where with p a real constant. 

Later on, we often write </>,^> for (j> u (x — cjt),ipu{x — u>i) when it is in a clear context. 
Substituting into the system (3) and demanding a solitary solution, we have 

H v (1, /x) — pH u (1, /x), 

and 


(15) - ujcj)" + (u> - 1)0 - H u (l, (i)^ 1 = 0. 

By (14), the solution of (15) is found to be 


(16) 


(p + 2)(w-l) 2 p Iuj 1 




Linearizing (3) around $ = (<f>, n </>) 4 , we have 

d 2 


where / 2X 2 = 


1 0 
0 1 


£ — ^ — UJ 'dx 2 _ ^ 2x2 

and \7 2 H U = 


-V 2 1A 


diuuijt* 1 VO di uv (<^>, VO 
H uv ((f>, VO H vv ((j),ip) 


Now our mission in this work is to 

(1) analyze £ and prove 

• it has exactly one simple negative eigenvalue A_; 

• the eigenvalue 0 of £ is simple; 

• the continuous spectrum of £ is bounded below by a positive number. 

(2) show d"(w) > 0 for d(w) = 0(<1> W ) +wf2(4> w ). 


2. Analysis on £ 


2.1. Eigenvalue and unitary transformation. We have 

vy2 TT _ f ddliuifillp) H UV iY ) l' l Y) A 


with 


and 

(17) 


_ ip f dd uu {\, p) dd U v( I5/1) \ 

V HuviX^) H vv { 1,/x) ) 

= V>oA4, 

V>o = ^(£ + l)(£ + 2)(w-l)sech 2 (|y^^a;), 


1 f dl uu {\, /x) H uv (l,/j,) A 

(p + l)iL u (l,/x) \ H uv (l 1 / 1 ) H vv (l, n) J 


We expect the existence of a unitary transformation, by which we can simplify 
A4 so that it is easier to analyze the spectrum, noting the unitary transformation 
does not change the spectrum of an operator. Since Ai is real symmetric, it suffices 
to find the eigenvalues of A4; if we can find two eigenvalues of A4, we find the 
unitary transformation matrix. 


Lemma 2.1. The matrix M. has two eigenvalues, which are Ai = 1,A2 = detAl. 
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Proof. From (17), we get 

det(Af - XI) = 


det (Ad) = det 


Huu — A(p + 1 )H U 


det (.Ad), 


//„ 


H„ 


Hw ~ X(p + 1 )H U 


We use the following result which shall be proved presently in Lemma 2.2 


Huu d~ pH uv (p T 1 )H U , 

H uv T pH vv — (p T1 )pH u . 

Multiplying the second column of Ad and adding to the first one, we have 

det (Ad — XT) = - - (1 ~ A) , det ( 1 „ V 

(p + 1)H U (1, p) \ p H vv A(p + 1)H U J 

Solving det(Ad — XI) = 0, we get Ai = 1, A 2 = det(Ad)/Ai, as the lemma states. 


□ 


Lemma 2.2. For H defined above, we have 


Huu (1; p) T pH uv (1, p) — (p T I)H U (1, p ), 

and 

Huv (1; p) A pH vv (1, p) — (p F l)pL7 1i (l, p). 

P+2 

Proof. Assume id(w, v) = CiU l v p+2 ~ l . The result can be easily verified by a few 

i =0 

differentiations and substitutions. □ 


Now we are able to form an orthogonal matrix O , by the normalized eigenvectors, 
such that 

°‘ M0 = ( J del^M) ) ' 

hence 

+ «(*!))' 
with <fo(z) = (P+ 1 )(P+2)(^-l) sech 2^^rT 2 ^ 

2.2. Spectrum analysis. We want to find out under what condition on Ad, oper¬ 
ator £ satisfies the spectrum conditions stated at the end of the previous section. 

Since O is an orthogonal matrix, it suffices to study C = O l CO. We want to 
prove 

Theorem 2.3. C is a self-adjoint operator on L 2 (IZ) x L 2 (1Z). It has a unique 
simple negative eigenvalue when det(Ad) < 0 the only other eigenvalue of 

C, which is simple; the continuous spectrum of C is bounded from below by a positive 
number. 


We shall complete the proof of Theorem 2.3 by showing a series of lemmas in the 
following. 

The self-adjointness of C is straightforward, since it is a compact perturbation 
of a selfadjoint operator [12]. Denote £ = ( ^ I with 
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First, observe that 


£2 = + ( “ 


1) — det(M)(j)o. 


AVb = Ai^i =► C ( ^ ) = Ar ( ^ ) , 

:* = A*^(»)=A 2 (;j. 

Namely, any possible negative (or 0) eigenvalue of C\ or £ 2 would make a (negative 
or 0) eigenvalue of C. 

We will show presently that C\ has a unique simple negative eigenvalue and a 
simple 0 eigenvalue, the only possibility for C to share the same property is that £ 2 
has strictly positive spectrum. 


2.2.1. Analysis on C\. Recall (fo = 1 ^ sech 2 (|hence 

r d 2 , ( ^ (p + l)(p + 2)(w-l) 2,P \u - 1 ^ 

£l = “"si + <" - 11 -2- sedl <2 V — x) 

is a compact perturbation of a self-adjoint operator and it is self-adjoint too. Also 
note this is exactly the linear operator we have by differentiating (15), where we 
got a byproduct, namely the eigenpair of C\: 

J eigenvalue : 0, 

) eigenfunction : <t>' {x). where </> is from (16). 

2 l -1 

Note that </> has shape of [sech (x)]p. Lemma 1.4 applies to (hence to C\ 
since u > 1). Therefore, we have shown 


Lemma 2.4. The self-adjoint operator C± has a unique negative simple eigenvalue 
and another simple eigenvalue 0. 


Moreover, we can compute the exact negative eigenvalue of C\. 
Let y = | then 


C\ — —i + ( w — 1) — 


p 2 (uj — 1) d 2 


4 u> 


p 2 {u} - 1) 


(_df_ 4_ 

V dy 2 v 2 


+ (w-l)-(/>o 

2 ) 


dy 2 

4_ _ 2{p+l)(p 

P' 


seclT 


p* 


It is easy to see that studying spectrum of C\ is equivalent to studying that of 


C-iv — —rrr H —T 


dy 

And it suffices to investigate 

•ClyO = — 


4 2(p+l)(p + 2) 


p* 


sech 2 (y). 


dy 2 


2(p + l)(p + 2) 


sech 2 (y), 


in which case the previous knowledge of Schrodinger operator comes into play. 

We know that C\ y has a unique negative eigenvalue. C\ v o might have more than 
one negative eigenvalue, but only the least one of C\ v q corresponds to the unique 
negative one of C\ y , with difference p-. So we focus on calculating the least negative 
eigenvalue of £i y o- 
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Comparing C\ v q with (5) and (7), we have 

2(p+l)(p + 2) 


a = 


a = s(s + 1) => s = 


1 | /I | 8 (p+l)(p+ 2 ) 

-t + i/IH - p + 2 


2 p 

By (9), the least negative eigenvalue corresponds to the largest e„, which is s 
when n= |V|. Namely we have the negative eigenpair: 

o 

Aiyo — —e„ — —s — —(1 + -) , 
ipi y o = -P s s (£ = tanh(y)) = sech s y. 

The eigenfunction comes from (10) and (11) with e n = s,£ = tanh(y). 

Now tracing back to L\ y and C\, we get their corresponding eigenvalues 

4 4 

A iy =Ai y0 + — = -(! + “). 

, _p 2 (w-l)^ -p(p + 4)(w-l) 

Al -- 4 - A iy -- 4 -■ 

2.2.2. Analysis on £ 2 . Our goal here is to find conditions under which £2 has 
strictly positive spectrum. Recall 


£2 = “ w ^2 + ( w - 1) - det(Ad)0o 

= -jAtizAl 1 + („_!)_ det(A<)0„ 


p 2 (w - 1 ) 


4w dy 2 

d 2 4 2(p + l)(p + 2) 


/ a 4 

V dv 2 v 2 


det(Ad)sech 2 (y)^ 


dy 2 p “ p- 

Lemma 2.5. TTie operator £2 has strictly positive spectrum if 

det(Ad) < ——-. 

p + 1 

Proof. It suffices to prove the same result for 

£ 2 „ = --+ 4 - 2(P + + 2) det(A4)sech 2 (?/). 


dy 2 p 2 


From lemma 1.3, a(C 2 V ) = ct c (£ 2 y) U cr p (£ 2 y) and a c {C 2 V ) = [4r,oo). So what 
we need to do here is to find what makes all the eigenvalues strictly positive. 
Similar to the previous subsection, for 


r r 4 d 2 2(p+l)(p+2) 2 

£ 2y o = £ 2y -n = -~r-o - 5 -det(A4)sech (y), 

dy z 


we have the least negative eigenvalues 

\ 2 j /0 


\ 2 
A 2 „n — — s , 


where 


-l + ^/l + ^ p+1 ^ p+2) det(Al) 
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And what we want is \ 2 V = — s 2 + > 0, which forces the following 


s* < — => s < - 

p z p 


^_ 1+ / 1 + 8 t£+l^+ 2 ) 4 

\ P P 

=> 1 + «E±WP + 2 l det(M) <( i + 1} 2 
=> 8 (p + l)g, + 2 ) det(A()< 8 ^ 2 ) 


det(TW) < 


1 


(P + 1)' 

We complete Theorem 2.3. □ 

Next, we will show Theorem 1.1 by finding out when d"(ui ) > 0 holds. 

3. Calculations for d"{uj) 

Theorem 3.1. We have the following for d"{ui): 

(1) If p < 4, then d"(u>) > 0. 

(2) If p > 4, then there exists a ui p > 1 such that d"{uf) > 0 for ui > u p and 
d"(co) < 0 for 1 < lo < ujp. 

Proof. We know the solitary wave solution = ((fuijifuY is a minimum point for 
the Lyapunov functional 

L(U) = G(U)+uSl(U). 

Hence 
Recall 
So 

Thus, 


£'($«) = ©'($*)+ wfi'($ u ) = 0 . 

d(pj) = 0($ u ) + wn($ w ). 


d\uj) =< L'($ w ), — 4> u > = n($ w ). 


d"(u) = — 


Using (fjijj for we shall investigate 

= \ f 4? + (t>ldx. 

A Jk 


Recall ft = 
we rewrite 


( p + 2 ) (^-i) S ech 2 (g- 


l -x) 


■ Let A - {P 2 + h!{Z)^ B = h/^ and 


cj) = Ap sech 2 {Bx) 


Z 1_ 2i 1 

cj) x = - ApB sinh(Bx)sechp + (Bx). 

P 

Sine if = pcf), we have 




1 + p 2 


(jr + (jr x dx 


m 
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i + p 2 


ApB / sech p(y)dy-\ — -pApB / sinh' {y)sec\\p + ~~(y)dy 


m 


m 


2 / id 2 id — 1 

= d i {to -1)?./-- + d 2 (ui - l) p \ -, 

V u> — 1 V u> 

where i?i > 0 , d 2 > 0 are constants. 

So 

d "(w) = 

= (oj — 1 )p~^oj~" s q{uj), 
where g(w) is a quadratic polynomial 




-($i + i9 2 )w 2 
P 


,01 2 ^ 
1 2 + p 


$ 2 

t ) w_ y 


We have trivially g(l) = 2^(i - ±) and 9 '(w) = |(0i +0 2 )w - ^ f )• 

So when p < 4, we have g(l) > 0 ,q'(oj) > 0, hence q(uj) > 0 and d"{w) > 0 for 
u> > 1. 

If p > 4, then g(l) < 0. For quadratic polynomial q(x) with positive leading 
coefficient |(i?i + $ 2 ), there is zero c p > 1 exists such that g(w) > 0 for w > oo p and 
q(w) < 0 for 1 < w < uj p . This concludes the lemma. □ 


Theorem 1.1 is concluded in view of Theorem 2.3 and Theorem 3.1. 


4. Examples 


Example 4.1. 

Ut + U x - u xxt + b\u p u x + b 2 {u p v) x + b 2 v p v x = 0, 

vt + v x - v xxt + biv p v x + b 2 u p u x + b 2 (uv p ) x = 0 . 

This is an example constructed in [15], with cq = l,ai = l,a 2 = 0. In this case, 


H{U) = b 1 - 


•i p + 2 


b\- 


i,p + 2 


'/P + 1 


-v + b 2 u- 


,,P+! 


(p + l)(p + 2) (p + l)(p + 2) p + 1“ ' "^~p+l' 

In [15], the author consider the proportional solitary wave with p = 1, finding the 
conditions for stability as 

(18) — b 2 {p + 2) < bi < —b 2 (p — 2), 


(19) 


— b 2 (p + 2 ) < &i < —b 2 


(p + 2)(p 2 — P — 4) 

p 2 + 3 p 


The method in this note gives a stronger result; the upper bound is relaxed. For the 
corresponding case where fi = 1 , using our method, one finds 


M = 


_1_ / H uu { l,p) H uv {l,p) \ 

(jp A H u (\, p) y H UV {1, p) Hw{ 1,m) / 

_1_ / h + pb 2 2 b 2 \ 

b\ + (p + 2 )b 2 \ 25 2 &i + pb 2 ) 


det(Al) < ——- => (p+ l)[(6i +pb 2 ) 2 - 46|] < (&i + (p + 2)5 2 ) 2 

p+ 1 

=> (p+ l)[( 6 i + (p- 2 ) 6 2 ][(&i + (p + 2)5 2 ] < (i>i + (p + 2)5 2 ) 2 

=> (p + 1 )[(&i + (p - 2 ) 6 2 ] < (&i + (p + 2 ) 6 2 ) 
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=> pbi < ~{p A l)(p - 3)6 2 
=> b, <— ijl+1 )( p — 3 > b 2 . 


To get the third inequality from bottom, the assumption H u = b\ + (p + 2)62 > 0 
is used, which gives a lower bound for bi, same as (18) and (19). To see that our 
result is stronger, first note that 62 > 0 holds true by comparing the first and third 
terms in inequality (18) (or (19)). Then, it is easy to verify that the upper bounds 
given here are larger than those given in (18) and (19). 


Example 4.2. Let q > 1 be an integer. Consider 


Ut 4“ u x u X xt A {v8v^ — 0, 

Vt ~ V x Vxxt A V^) x = 0. 

This is an example considered in [9] , where the result (Lemma 3.1) requires uj — 
1 + H uv — H u 2 > 0, which is obviously a much stronger assumption than ours. This 
is exactly the point: by taking advantage of the accurate point spectrum information, 
we weaken conditions like this. 

By our result 


H(u , v) = - 

Existence of solitary wave solution forces p 
We have here 


g + 1 

= 1. This case corresponds to p 


2q—l. 


KA _ _1_ f H uu [\,p) H uv (l, p) \ 

{p + p) \ H uv (l, p) H vv (l, p) J 

J_ ( q(q- 1) q(q + 1) A 
2 q 2 V g(g +!) q(q- i) ) ' 

Hence 

det(TW) = -- < 0 < —j—■ = - 1 -. 

q p+ 1 2 q 

The condition is satisfied, so Theorem 1.1 applies. 


Example 4.3. Let q > 1 be an integer. Let A, B ^ 0, consider the decoupled system 

Ut A u x u xx t + ( Au p + ) x = 0, 

Vt A v x Vxxt A (Bv pJr ) x — 0 . 


In this case 


H{u,v ) = 


AuP+ 2 + Bv p+2 


q A 2 

H v (l,p) = pH u (l,p)^Bp p+1 = 
=> p = 0, or Bp p = A. 


Ap 


A A _ _1_ ( H uu {\, p) H uv (l, p) \ 

(pAl)H u (\,p) \ H uv (\,p) H vv (l,p) J 

J_ _ ( A(p+ 1) 0 \ 

(pA 1)A V 0 B(p+l)p p J ' 

Hence 

det (M) = 
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If p = 0, then det(.A4) = 0 < Theorem 1.1 applies. 

If p ^ 0, then det(.A4) = 1 > 

This is a slight generalization of example in [5]. It tells that, the proportional 
solitary wave solution of a decoupled BBM system is generally not stable, unless one 
component is 0. Indeed, if the solitary wave is perturbed in such a way that 

the initial data are exactly a new pair of solitary waves (<j > Ul , i/; W2 ) that are close to 
(4>ui, ifui), but ^ W 2 , then the solution (new pair of solitary waves) will slip away 
form ,ip u ) as time goes on, hence not stable according to the definition we are 
using here. This seems due to the proportional requirements. This conclusion more 
or less goes against intuition. The fact that each component is stable while the 
combination is unstable actually brings up interesting questions: does there exist 
coupled non-proportional solitary wave solutions and if so, are they stable? Is there 
a different definition of stability that can accommodate the missing stability that 
solely caused by the phase difference? 

In the study of vector solitons in birefringent nonlinear optical fibers, the nonlin¬ 
ear Schrodinger equations have been proved numerically to have non-proportional 
solitary wave solution [19]. Next, the authors hope to investigate the similar subject 
in other dispersive equation system such as coupled BBM model. 

The final example is an analogy to the one in [19] for a coupled nonlinear 
Schordinger equations. 


Example 4.4. 

Ut + U x - u xxt + {(u 2 + /3v 2 )u) x = 0, 
v t + v x - v xxt + {{v 2 + 0u 2 )v) x = 0. 

In this case p = 2 and 

TT . u 4 0u 2 v 2 v 4 

H (u,v) = — H — + —. 

H v { 1, p) = pH u { 1, p)^(0- l)/x(l - p 2 ) = 0. 


If /3 = 1 , then det(.A4) = | = 

If 0 ^ 1, p = 0, then det(.A4) = f < =>■ /3 < 1. 

If 0 ^ 1, p ^ 0, then p 2 = 1, det(.A4) < ^ => 0 2 > 1. 
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